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This paper presents a second-order statistical analysis for the cascade of 
a bandpass limiter, and ideal phase detector and a video filter. This cascade 
forms an important subsystem in the mathematical model of some coherent 
communication systems where information is transmitted by phase or 
frequency modulation of the carrier. We derive the autocorrelation function 
R{t\ , ta) of the video filter response when the bandpass limiter input is a 
fixed amplitude-phase modulated carrier plus stationary gaussian noise. 
The video filter response is wide sense stationary for some nontrivial cases; 
these include biphase, single tone, and stationary gaussian noise phase 
modulation. For these cases, we obtain the video filter output average power 
spectrum as the Fourier transform of R(t) for all values of the limiter input 
signal-to-noise power ratio. An application of the results of this paper is 
the performance of a FM-PM demodulator for a set of parameters charac- 
teristic of one mode of operation of the Apollo Unified S-Band communica- 
tions system. We present the performance as a family of curves of subcarrier 
channel output signal-to-noise power ratio as functions of the limiter input 
signal-to-noise ratio where subcarrier phase modulation index is a param- 
eter. The approach is similar to the analysis by Davenport of the signal- 
to-noise ratio transfer characteristic of an isolated bandpass limiter. 

I. INTRODUCTION 

In some coherent communication systems, such as the Apollo Uni- 
fied S-band system, 1 where information is transmitted by phase 
modulating a carrier, bandpass limiters 2 are used in the IF channels 
preceding the coherent demodulators. Ideally the bandpass limiter 
removes any amplitude modulation that might exist before the signal 
is demodulated. 

233 



234 



THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1969 



Figure 1 shows a typical coherent phase demodulator used in such 
a system. This demodulator consists of a multiplication operation (a 
phase detector) with post-video filtering. The phase modulated signal 
is multiplied by a coherent carrier reference to yield a video signal 
containing the desired information. The signal into the limiter is 
usually accompanied by noise that is frequently assumed to be addi- 
tive and gaussian. The presence of the noise affects the performance 
of the demodulator in a very complicated way because of the non- 
linearity of the limiter. Thus it is difficult to evaluate the corruptive 
effect of the noise on the demodulated information. 

One criterion of performance at points in a communication system 
is the signal-to-noise power ratio (S/N). For the cascade in Fig. 1, a 
problem of interest to the systems engineer is the video filter output 
S/N as a function of the input S/N to the limiter when the input noise 
is additive, stationary, and gaussian. The relationship is known between 
input and output S/N for an ideal bandpass limiter where the input is 
the sum of a stationary gaussian noise and a signal P(t) cos (w e t + <j>) 
(see Ref. 2). For the analysis there, P(t) is a random process and is 
slowly varying compared with cos u e t. The carrier phase <f> is a random 
variable independent of P(t) with a uniform distribution over [0, Shr]. 

It is not possible to apply the known S/N transfer characteristic of 
the ideal bandpass limiter found in Ref. 2 directly to obtain the S/N 
transfer characteristic for the bandpass limiter-phase detector-video 
filter cascade. A knowledge of the form of the signal and the noise 
out of the bandpass limiter, and not just the S/N of this output, is 
necessary to determine the effect of the phase detector on the bandpass 
limiter response. 

To obtain the cascade S/N transfer characteristic we apply the 
mathematical tools used in Ref. 2. The form of the signal assumed in 
the analysis of the cascade is s(t) = P cos [co e + 0(t) + <t>] where P is 



P (t) COS [ci» c t + B (t) + <p] + NOISE 
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Fig. 1 — A coherent phase demodulator with IF bandpass limiting in the 
presence of additive noise. 



DETECTOR- FILTER CASCADE 235 

a positive constant, 6(t) is phase modulation that is slowly varying 
compared with cos u e t, and </> is a random variable representing the 
arbitrary initial phase of the signal carrier. The probability density 
function of <j> is assumed to be uniform in the interval [0, 2ir]. The 
noise input to the bandpass limiter is assumed to be additive, stationary, 
and gaussian with zero mean and power spectral density N. The input 
noise, the modulation 6{t), and the carrier phase <f> are assumed to be 
jointly statistically independent. For the following analysis, the limiter 
is assumed to be ideal with limit level I. The transfer function of an 
ideal limiter is denned by 

[ + Z, x > 

y = l( x ) =\ 0, x = (1) 

[-1, x < 0. 

A coherent carrier reference sin (a e t + 4>) is assumed to be available 
for the demodulator where </> is the phase of the carrier. 

II. THE SECOND ORDER STATISTICAL ANALYSIS 

2.1 A Cascade Model when s(t) is Narrow Band Limited 

In order to obtain a S/N transfer characteristic for Fig. 1, the auto- 
correlation function of z(t) is derived. When R^k , t 2 ) = R,(t) the aver- 
age power spectrum of z(t) is denned by the Fourier transform of R.(t) 
and the S/N transfer characteristic can be found. An analysis of the 
autocorrelation function of z(t) does not seem possible for general s(t). 
However, if the signal s(t) is a narrow band-limited process such that 
the bandpass niters are narrow compared with the carrier frequency io e , 
the response z(t) should be the same with or without the post bandpass 
filter that precedes the phase detector. The response of the nonlinearity 
l(x) to an input x(t) = s(t) + n(t) that is narrow band-limited about 
±a} e is a family of terms narrow band-limited about the frequencies 
±nw e where n = 0, 1, 2, 3, • • • (see equation 13-53, section 13-1 of Ref. 3). 
Any narrow band-limited input to the phase detector that is not about 
±w e will generate a phase detector response above the cutoff frequency 
assumed for the video filter. For a narrow band-limited x(t) the auto- 
correlation function of z(t) is obtained from the analysis of Fig. 2. 

2.2 The Derivation of the Autocorrelation Function of z(t) 

Assume that the input x(t) is narrow band-limited such that Figs. 1 
and 2 yield equivalent z(t). The autocorrelation function R.{t x , t 2 ) is 
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Fig. 2 — The narrow band equivalent receiver for the derivation of R,(t\, t t ). 

obtained by first deriving R a (t x , t 2 ) from the model in Fig. 2. Since z 
and w are related by the linear video filter, R^U , t 2 ) follows directly 
from R W {U , t 2 ). 

The Laplace transform solution of a zero memory nonlinearity with 
stochastic excitation is used to derive R w (ti , t 2 ) (see Chapter 13 of Ref. 3). 
The limiter characteristic is 



m 



= 2~'\ J '*^ 6XP ^ dl ° ~*~ J ^"^ GXP ^ do> 



(2) 



where 



and 



r +0 ° i 

/+(w) = / l(x) exp (—cox) da; = - , for Re [co] > 
Jo w 

r° I 

/_(co) = / Z(x) exp (—cox) dx = - , for Re [co] < 0. 

J -co w 

The variable co = w + ;w is complex with Re[co] = u. The contours 
C+ and C_ are taken parallel to the v axis in the co plane with Re [co] > 
for C+ and Re [co] < for C_ . For convenience l(x) is written symbolicly 
as 



Kx) = g— . J /(co) exp (oca) da 



(3) 



where equation (3) means the same as equation (2) when C t and C_ are 
not the same contours. 

Since w(t) = sin (co c f + <f>) -l[x(t)], the autocorrelation function of 
w(t) is 

R(U , t 2 ) = [y^jj J 1M J f(co 2 )#{sin («.(, + 0)-exp (co,s, + co,n,) 

•sin (co e £ 2 + <£)-exp (co2S 2 + co 2 n 2 )} cfco, dco 2 (4) 
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where s,- = s(U) and n, = n(t { ), i = 1, 2. The order of complex integra- 
tion and the expectation operation have been interchanged to get equa- 
tion (4). For the assumed statistical independence of n(t), 0(0, and <f>, 
the expected value in equation (4) factors into 

22 {sin (uJi + <£) • exp (wtfO -sin (u c t 2 + <f>) • exp (w 2 s 2 ) } 

• exp $[<r 2 «? + 2R n (T)u l w 2 + <r 8 «a] (5) 

where t = t 2 — ^ . The form for the cross correlation function 
.Ejexp (wxTix) exp (0^2)} where n(t) is stationary gaussian noise has been 
used in equation (5) (see pp. 476-477 of Ref. 4). 

For the case where s(t) is narrow band-limited with respect to w e , 
the filter in Fig. 2 is a narrow bandpass filter, and is assumed to be 
symmetrical about w e . Then n(t) can be written as (see pp. 373-374 of 
Ref. 4) 

n(t) = x c cos io c t — x. sin u) e t 

where x e and x. are statistically independent stationary gaussian random 
processes, and 

R„(t) = R,(t) cos u c t (6) 

where R,(t) = R Xc (t) = #».(*■)• For a narrow bandpass IF filter, the 
transform of R r (r) is lowpass with a narrow bandwidth compared to u e . 
With the substitution of 

ti = t, t 2 = t + t, 

4>* = $ + u c t, 



.,. _ e xp (#*) - exp (- # >*) 



sin«* = 2 . 

and 

+ 00 

exp [ J B n (T)co 1 o) 2 ] = 2 /«(uiwjfi») exp (jmu c T) (7) 

B1--O0 

(see Article 1, Chapter 3 of Ref. 5), equation (5) becomes 

+ 00 

(— 1) 2 JmttW^) exp (jwiw e T)-I?{[exp (jco c t + j20*) 

m — — oo 

+ exp (— jw e T — j20*) — exp (ju e r) — exp ( — jw c t)] 

• exp [co,P cos (0, + <(>*) + co 2 P cos (d 2 + <*>* + w.r)] } . (8) 
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Since exp 0'w e r) and exp [w 2 P cos (0 2 + <£* + w c t)] are periodic in u c t 
with the period 2ir, the function 

exp (jmu e T)E{ } (9) 

in equation (8) is periodic in w t. Since R„(t) transforms to a narrow 
band-limited lowpass spectrum, the autocorrelation function of z(t) 
corresponds to the dc component of the Fourier expansion of equation 
(9). With the substitution of S = w c t, the dc component of equation 
(9) is 

J — exp(jmS)-E[ } 

= E E IMP)-I k (c 2 P)-E B < f. 

■[^♦{exp [/(w + 1+ fc)5 + i(2 + r + fc)0* + ;M + kd 2 )]} 

+ ^.{exp [;(m - 1 + k)S + i(-2 + r + fc)<£* + j(r0, + k6 2 )]\ 

- E,.{exp \j(m + 1 + k)8 + j(r + fc)0* + j(rtf, + fc0 2 )]} 

- £,.{exp [i(m - 1 + k)8 + i(r + &)<*»* + ;'(r*. + fcfc)]} j|- (10) 

Since <t>* = u e t + <£, <f>* has a uniformly distributed probability density 
function on [0, 2ir]. The averages in equation (10) with respect to 5 
and 4>* follow. For example, the first average with respect to 5 and <f>* 
is zero ifco + fc + 1 ^0or^ + r + 2?^0, and when k = — 1 — m 
and r = — 2 — k = m — 1 the double average is exp [(m — 1)0, — 
(m + 1)0 2 ]- Equation (8) reduces to 

(-i) E UwM] lL a -\ F '/Lr^D-Efexp Km - 1)^! - j(m + 1)0J| 

m=— oo L 

+ I { m VPl-Z P -\yE{exp [j(m +1)6,- j(m - 1)0,]} 

- liVi P) I l -?J?»E{exp [j(m + 1)6 X - j(m + 1)0 2 ]\ 

- iL^riLT^Eiexp [j(m - 1)0, - j(m - 1)0J}]. (11) 

The terms in equation (11) for positive and negative m can be combined 
by noting that I- m (x) = I m (x). With the substitutions 

uwm = z 2m+ i i ql ; ( 2 m + - ff + 1} d2) 
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and 

h m , k = -^ ffoyf exp iff)l m {*P) dco, (13) 

the autocorrelation function of z(t) is 

BXtx , Q - \ E E 2 9 « + -g!(g + m)! 

•[/t 2 ra+ i,2, + X(m + 1, m + 1, *i , fc) 
+ J&-i.**«R*(>n - 1, » - 1, t , y 

- A H +i.»t+Jk*-i.u+-Bf(ifi + !j m — !i *i » *») 

- fc. + i.i« + ^-- i..t+-Bf(« - 1, m + 1, «, , Q] (14) 
where 

fl, m = 
1.2, m > 
and R e (A, B, *, , fe) = #{cos [A0(i,) - B6(t 2 )}} for any integers A and £. 

III. THE CLOSED FORM SOLUTION FOR h m _ k 

The autocorrelation function of z(t) given in equation (14) contains 
the constants h m , k where m + k are odd integers. For the ideal limiter 
characteristic of equation (1), there are closed form solutions for these 
parameters. Since /+(«) = l/w for Re [w] > and /_(«) = l/w for 
Re [co] < 0, equation (13) becomes 

*-■» = 2^- / c /a,t " 1/ "' (coP) exp (t") rfw 

+ 27-/ c /w *" ,/ -( wp ) exp (t~) rf " (15) 

where C_ is the contour (— e — j», — e + j«) and C + is the contour 
(+e — j oo, +e + joo). By the change of variable w = jaj and the sub- 
stitution of I m (z) = (j)~ m J m (jz), analytic continuation can be applied 
for m ^ and k ^ to give 

K. h = l - Cf)**- 1 /" x'*-" J„(xP) exp [-=1^] ^- (16) 

When m + fc is even, the integrand of equation (16) is odd and h m , k = 0. 
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When m -f- k is odd, the integrand of equation (16) is even and 

( m + k\(py 2 

21 k+m ., l \ 2 l\2a 2 ) ( m + k . m , t . z^l) rm 

h m , t = - (j) r-j k iF l [— T - ,m + 1,-gr) (17) 

2T(m + l)[ W \ 
where a solution has been used for the integral 

p x k - l J m {xP) exp z= y^J dx (18) 

in terms of the confluent hypergeometric function \F\ifit\ /3; —x) (see 
equation A. 1.49, p. 1079 of Ref.6). For the case when m and k are non- 
negative integers iF^m + k/2; m + 1; —x) can be expressed in closed 
form in terms of first and second kind modified Bessel functions. A 
list of these expressions is given by Middleton (see equation A. 1.31, 
section A 1.2 of Ref. 6). A collection of h m>k in closed form for low order 
indices is given in Table I. For Table I, x = P 2 /2a 2 is the input signal- 
to-noise power ratio into the limiter in Fig. (2). 

Any of the h m _ k in equation (14) can be found in closed form from 
Table I by using the recurrence relations 

_ h 2(ro + 1) , , 4(m + l)ro , , 1Q v 

lt m +2,k — n m.k — p ">m-l,k-l ~T p2 llm.k-2 > V ly / 

P (fc - w - 2) 

'•'m + i.k+i — 2 "-m.k 2 »*m-l.fc-l 

<T 0" 

+ 2(t- ? -a)» > _ (20) 

and 

(m - fe) , P 2 , (ro - k) . . 

h m ,k + 2 = 2 hm.k H 4"m-2.ft 4 -T«m-1.*-1 . (.21; 

CT (T (T 

Equation (19) is derived from equation (16) by using the Bessel 
function identity 

J m+2 (xP) = 2{m + 1} J m+1 (xP) - J m (xP). (22) 

Equation (20) is derived through a by-parts integration of equation 
(16) and the application of equation (19). Equation (21) is derived 
through by-parts integration of equation (16). In the development of 
equations (19), (20) and (21), the integral in equation (16) is re- 
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Table I — Closed Form Solutions of Some h m , k 
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stricted to the half interval [0, oo) which is possible since the inte- 
grand of equation (16) is even when m+fc is odd. 

IV. THE AVERAGE POWER SPECTRUM OF z(t) 

The autocorrelation function of z{t) given in equation (14) becomes 
time independent such that z(t) has the average power spectrum S,(u>) = 
F[R m (t)] when R {A, B, t x , t 2 ) = Rg(A, B, t) for integers A and B. 
There are some important cases of 6{t) for which R e is time independent. 

If 8 is a biphase modulation with 6{t) = ± | | that has a zero mean 
and autocorrelation function (see equation 9-42, section 9-2 of Ref. 4) 

\\ 6 | 2 (l - Jfi) , for | r I ^ T 

flt(r) = 1 (23) 

0, for | t | > T 

then 

Rg{A, B, t t , t 2 ) 

= cos A | 6 | • cos B | 6 | + sin A \ $ \ -sin B \ 6 \ -r t (r) (24) 

where r e (r) = R (j)/\ 6 | 2 is the normalized autocorrelation function 
of d(t). Then R e is a function of t = t 2 — t\ . 

For a single tone modulation given by 8(t) = mj sin (co,£ + £) where 
£ is a random variable with a uniform probability density function 
on [0, 2r]| a simple Bessel series expansion gives 

CO 

R B (A, B, t x , U) = 2 ^nJ2n{Am^J 2n {Bm^ cos (27kojT) 

n-0 

+ Z *nJ2«-iUm l )J in - l (fim 1 ) cos [(2n - 1)«,t] (25) 
n-i 

where 

|1, n = 

l2, n > 0. 

For the single tone modulation, R $ depends only on the time difference r. 
If 6(t) is the sum of tones 

0(0 = £m p sin(co p « + £ p ) (26) 

p-i 

where £ p , p = 1, • • • , N, are independent random variables with 



(27) 



(28) 
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uniform probability density functions on [0, 2v], Re is again independent 
of time. 

If 0(0 is a stationary gaussian process with zero mean, variance o> 
and autocorrelation function K (t), then Re(A, B, t x , t 2 ) = Re(A, B, t). 
The second-order characteristic function for the stationary gaussian 
process is denned as (see equation 112, Chapter 7 of Ref. 4) 

*•(«! , « 9 J r) = .E(exp |j[<M(« + t) + w a »(«)]}) 

= exp [-K»(0)(w? + <&) - £ # (r)«,wj. 
Then 

fl,(A, B,t x ,t 2 ) = Real Part E{ exp tf-M - i^0 2 )} 

- exp [-| (A 2 + B 2 )]-exp [ABK B (r)] 

= R e (A,B,r). 

The validity of equation (14) depends on the narrow band-limited 
assumption for the modulated signal s{t) at the carrier frequency 
w . For s(t) to be narrowband limited, the parameter values that the 
modulation functions can have are restricted. 

V. AN APPLICATION OF THE it!- RESULTS TO THE PERFORMANCE OF A SUB- 
CARRIER CHANNEL 

A modulation technique sometimes used for communication is FM- 
PM where the carrier is phase modulated by a subcarrier that is in 
turn frequency modulated by the information waveform. The FM- 
PM signal is of the form 

8(0 = P cos [a c t + + m x sin [uj + £ + \(t))\ (29) 

where P, a> e , w, and m, are constants, <f> and £ are independent random 
variables usually assumed to have uniform probability density func- 
tions over [0, 2tt], and \(t) is the integral of the information waveform. 
In a typical application, w e » w, and \(t) is slowly varying compared 
with cos oji<. With these restrictions the information X(0 can be re- 
covered from s(t) with the receiver shown in Fig. 3. 

The purpose of the bandpass limiter is to remove the effect of varia- 
tions that might occur in P. For the ideal case where s(t) is not per- 
turbed by noise, the subcarrier filter input 2(0 is 



244 



THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1960 



S(t),N 


BAND- 
PASS 

FILTER 


t 
1 


Mi) 


T 
i 


BAND- 
PASS 
FILTER 


r®^ 


VIDEO 
FILTER 


t 

1 


SUB- 
CARRIER 
BAND- 
PASS 
FILTER 




SUB- 
CARRIER 
DEMOD- 
ULATOR 






x(t) 


yet) 




zCt) 







REFERENCE 
SIN (Ctf c t+tfO 



Fig. 3 — FM-PM receiver with ideal bandpass limiter. 



21 . 



z(j) = sin {mj sin [uj -j- \(t) + £][ 



41 



(30) 



= E J 2 n- i (m i ) sin \{2n - l)[co,* + \(t) + £)}. 



If X(0 is slowly varying compared with cos co,£, the information can be 
recovered with a subcarrier filter that passes only the first component 
of the sum in equation (30). For the noiseless case the subcarrier filter 
response is then 



-- J l (m 1 ) sin [w,t + \(fl + £]. 



(31) 



After additional processing in a subcarrier demodulator, X(t) is obtained 
from equation (31). One criterion of performance of the receiver is the 
S/N out of the subcarrier filter as a function of the limiter input S/N, 
x = P 2 /2(t 2 . Since \{t) varies slowly compared with cos uit, the output 
S/N for the subcarrier filter is determined with sufficient accuracy by 
setting \(t) = 0. If \(t) = 0, the subcarrier output S/N follows directly 
from equations (14) and (25). Substitution of equation (25) into equa- 
tion (14) gives the power spectrum 



S.(«) 



S 2h\ Q 2 / 2 2 „-i(m 1 )-F[cos (2n - 1)«,t] 



■ in («,! + {) 

+ {l)\*K - <Th 21 J (2 mi )] 2 -F[r v (T)] 

+ {hWhl, J2 Jl(2m 1 )-F[r,(T)-oo8 (nw.r)] 

n-l 

+ (4)<r 4 A? 2 2 J r 2 2 n - 1 (m 1 )-F[r^T)-cos(2n - l)« lT ] 

n-l 
eo 

+ (tV) £ *[o**i»«/.(»0 - ff 2 Asa^«(3m 1 )] 2 -F[r*(T)-cosnco 1 T] 
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+ (&)[A» - ^ 23 ^o(2m 1 )] 2 -nr'(r)] 
+ (A)*X E JlVrnJ-FirKr)- COS m, lT ] 

n-l 

+ (Fir) Ee„[<r 3 / l2 3^(2m 1 ) - a% i3 J ^m,)] 2 ■ F[t\{t) ■ cos ruo.r] 

n-0 

+ (&)<rX S J 2 2 n -,(m 1 )-nr:(r)-cos(2n - l)« x r] 

n-l 

+ (rkr) E eJiaXMrnJ - <7 4 / l34 J n (3m 1 )] 2 -F[rt(T)-cosrto 1 T] 

n-0 

+ Oriir) Z «»[>W«(3»ii) - crX^Sm.^'-n^C^-cosTia)^] (32) 

n-0 

where r„ = R,/R,(0) = RJa 2 . The approximation, equation (32), 
neglects all the terms of equation (14) containing the factor Rl a+m 
where 2q + m > 4. The terms in equation (32) are the significant terms 
of S t (u>) for the single tone modulation. The spectrum in equation (32) 
is the weighted sum of terms of the form 

F[r n t (T) cos mo),r] = ^ *TC(t)] * F[cos ma> lT ] (33) 

where * is the convolution operation. Since F[cos mcoir] is a pair of 
impulses of weight -k at ±mwi , 

F\f v {j) COS wza>,T] = £[£».„(" + wa>i) + £„,„(w - Tnco,)] (34) 
where 

S,» = F[r:(r)]. 

The first term in the spectrum of equation (32) is the signal content 
of z(t). All other terms of equation (32) correspond to noise alone or 
a combination of signal and noise. All terms of equation (32) except 
the first term are usually combined to give the interference (noise) 
spectrum at the output of the video filter. 

A computation was made for the subcarrier filter output S/N as a 
function of the input S/N x. The following conditions are assumed for 
the computation. 

(i) The power spectrum of the input gaussian noise to the cascade 
in Fig. 1 is uniform over the bandwidth of the prelimiter bandpass 
filter. 
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(u) The prelimiter bandpass filter is assumed to have a gaussian 
transfer function such that 



[=7*1- 



r.(r) = exp \-^\- (35) 



[Hi) The subcarrier amplitude transfer function is 



Hfr) | = 



. Ac* . , , Aco 

lj W, — < I CO I < CO! + — 



(36) 



.0, all other u, 

where Aco « co . Also, co = 12.566 X 10 6 and w, = 6.434 X 10 fl are 
assumed. Substitution of equation (35) into equation (34) gives 

(37) 



ft-M-ra«-^-*[-s(yv 



From condition Hi, the noise spectrum in the passband of the subcarrier 
filter is approximately constant when co = co x . The signal and noise 
powers out of the subcarrier filter follow from S c (^i)- The signal power 
is 2h\ J\{m 1 ); the noise power is 

Y f [S.(w) - 2ti J 2 i(md-n™so> 1 T)]-\H(jcS) [ 2 tfco S [B'.M]-2 A/ 

where A/ is the width of the subcarrier filter and where S' t is equation 
(32) with the first term omitted. The function 

s ^ = SSfr (38) 

was computed for x between 0.01 and 100 with m^ as a parameter. The 
results of the computation are shown in Fig. (4) . For a given m, and x, 
the output S/N for the subcarrier filter is a:/2A/-*S(m 1 ). 

VI. SUMMARY 

A general, second order statistical analysis is presented for the cascade 
of a narrow bandpass limiter, an ideal phase detector, and a video filter. 
In this analysis, the input to the limiter is assumed to be the sum of a 
stationary gaussian noise and a fixed amplitude phase modulated sine 
wave. The autocorrelation function of the cascade response is obtained 
as a function of the signal-to-noise ratio x at the limiter input, the nor- 
malized autocorrelation function of the lowpass equivalent for the 
limiter input noise r,(r), and the phase modulation 6{t). 

The cascade response z{t) has the autocorrelation function R t {U > ^2) 
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Fig. 4 — The unit bandwidth subcarrier filter output, S/N normalized by x 
where 10" 2 ^ x ^ 10 2 and Mi is a parameter. 

that can be time dependent. However, for some important cases of 0(0, 
R,(ti > h) = R*(t)> and the cascade response has the average power 
spectrum S,(co) = F[R,(t)] where F is the Fourier transform operation 
with respect to t. The cases of 0(0 considered that yield R t (r) are the 
random biphase waveform = ± | |, the single tone 0(0 = w x sin (co^+S), 
and the stationary gaussian process with autocorrelation function K 9 (t). 

The dependence of R t (tx , t 2 ) on the limiter input S/N appears in the 
h parameters. These parameters can be obtained in closed form as func- 
tions of the modified Bessel functions h(x/2) and Ii(x/2). The lower 
order h parameters encountered in the first few terms of the series for 
R, are found, and recurrence relations are derived through which 
higher order h parameters can be derived easily. 

For the modulation types that make R, a function of t alone, the 
power spectrum S,(co) is known for all values of the limiter input S/N x. 
Then the S/N can be derived in any frequency band at the output of 
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the video filter in Fig. 1 as a function of any S/N into the limiter. 
The performance of a subcarrier channel was considered where 0(0 = 
m, sin [wit + X(0 + £]. The subcarrier was assumed to be phase mod- 
ulated by a narrowband low pass process \(t). The S/N at the output 
of the subcarrier filter was obtained by computation of the approxima- 
tion of equation (32). For this example, a gaussian prelimiter bandpass 
filter was assumed. For this filter shape, r"(r) and its transform S v . n (u) 
are gaussian for all integers n. Some representative parameters from 
the Apollo unified S-band communication system 1 were assumed. These 
were 

(i) A prelimiter noise equivalent bandwidth of 4 MHz. 
(ii) A subcarrier frequency of 1.024 MHz. 
(Hi) A subcarrier noise equivalent bandwidth of 0.2 MHz. 
(iv) An input S/N range of 0.01 ^ x g 100. 

(v) A set of modulation indices m, = (0.2)fc, k = 2, 3, 4, 5, 6, 7, 8, 
9, 10. 

The results are given in Fig. 4. The differential between subcarrier filter 
output S/N at low and high values of x is a monotonically increasing 
function of m x for 0.4 ^ m x ^ 2.0. The shapes of the curves are similar 
to that of the (S/N) /(S/N), curve obtained by Davenport. 2 
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